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Abstract 

We present an analytic calculation of electron and tau C(a 3 ) loop effects 
on the muon anomalous magnetic moment. Computation of such three-loop 
diagrams with three mass scales is possible using asymptotic and eikonal ex- 
pansions. An evaluation of a new type of eikonal integrals is presented in 
some detail. 



1 Introduction 



The new measurement of the muon anomalous magnetic moment, = (g ll — 2)/2, in 
the experiment E821 in Brookhaven has motivated many recent theoretical studies. 
In the Standard Model, a M receives contributions from electromagnetic and weak 
interactions, as well as from loop effects involving hardrons. All three types of effects 
have been studied recently. QED contributions are known to the four-loop level, and 
even some five-loop diagrams have been evaluated. The next largest contribution 
is due to hadronic loops and is the most difficult one to evaluate. There has been 
recently significant progress in both evaluation of the light-by-light diagrams [1] 
and hadronic vacuum polarization effects [2]. Electroweak two-loop effects are also 
known [3]. 

The present paper is devoted to the only 3-loop QED contribution to a M that has 
not been evaluated analytically so far: from a diagram with electron and r lepton 
loop insertions in the photon propagator, shown in Fig. 1. Because of the three 
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Figure 1: Three-loop contribution to muon anomalous magnetic moment with the elec- 
tron and tau loops inserted in the photon propagator. 

mass scales present in this diagram, m e ^ T1 it cannot at present be evaluated in a 
closed form. However, we present an approach based on asymptotic and eikonal 
expansions, which takes advantage of the wide separations between those scales and 
permits an evaluation with an arbitrary accuracy. 

We begin with a brief summary of the present knowledge of QED contributions 
to the electron and muon anomalous magnetic moments, a efM (for more details see 
e.g. [4]). 

(a) Electron. To match the present experimental precision, one needs four terms 
in the expansion of a e in the fine structure constant a, 

« e =E^(-y +■■■ (1) 

n=l W 

where ellipses indicate contributions of loops containing the heavy leptons \x and r; 
A± t 2 have been known since the early years of QED [5]. An analytical evaluation 
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of A 3 required the efforts of many groups and took almost 40 years; it has been 
completed only recently [6]. A 4 is known only numerically. 

The perturbative series for a e is very well behaved, and, together with the most 
recent experimental values for electron and positron [7], allows one to deduce a very 
precise value of the fine structure constant [8, 4], 

a' 1 = 137.03599959(38)(13). (2) 

(b) Muon. The value of a found from electron g — 2, eq. (2), can be applied 
to compute the QED contribution to the muon anomalous magnetic moment a M . 
Because of the presence of electron loops, higher-order QED contributions to are 
enhanced with respect to a e . At present five terms of the expansion in a are needed: 

with 

Ci = A x = 0.5, 

C 2 = A 2 + a^mjm^) + a 2 (m ll /m T ) = 0.765 857 388(44), 

C 3 = A 3 + CT\e) + C 3 77 (r) + C 3 vac - poL (e) + C 3 vac - poL (r) + C 3 vac ' poL (e, r) 

= 24.050 509(2), 

C 4 = A 4 + 127.55(41) = 126.04(41), 

C 5 = 930(170), (4) 

where ai )2 describe contributions of two-loop diagrams with electron and tau loops, 
respectively [9, 10]. We will discuss them in detail later on. 

In C 3 we have contributions from light-by-light scattering diagrams with e and 
r loops [11], and vacuum polarization diagrams with either e, or r [12], or both 
types of loops. An analytical evaluation of this latter contribution of mixed e — r 
diagram, C™ ' pol '(e,r), is the main purpose of this paper and will be presented 
below. Numerically it can be evaluated using the kernel from [13]: 

C vac. poi.^ r ) = 0.0005276(2). (5) 

For C4 one uses the difference between the muon and electron coefficients found in 
[14] . For C 5 only a numerical estimate of the presumably dominant contributions is 
known [15]. 

The present estimate of the total QED contribution to a M is [4] 

a QED = 116584705.6(2.9) x 10" 11 . (6) 
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2 Two-loop diagrams with electron and tau loops 



Before calculating the contribution of the three- loop diagram of Fig. 1, we would 
like to discuss a method of evaluating the two-loop diagrams shown if Fig. 2 (we 
discuss only the Pauli formfactor, relevant for a M ). It will serve us as an example to 
illustrate the main points of our calculational techniques. Of course, full analytical 
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Figure 2: Two-loop contributions of lepton loops in the photon propagator. 

results for these two-loop contributions are known [9, 10]. It is instructive to note 
how they were obtained. First, a closed analytical formula was found, with full 
functional dependence on the ratio of masses m e /m M or m M /m r [10]. That formula, 
containing dilogarithms, was found awkward to use, because of cancellations and 
difficulties in estimating the accuracy in the numerical evaluation. In [9] expansions 
of the exact result in powers of small mass ratios were given. Such expansions avoid 
evaluation of special functions and their accuracy can be precisely assessed. 

Here we demonstrate how such expansions can be obtained without knowledge of the 
exact result. We first consider the case of the r loop insertion. In this case a well 
known method of heavy mass expansion is applicable (for a review see [16]). Let us 
denote the loop momenta by p T and p M , for the momentum flowing inside the r loop 
and for that in the virtual photon respectively. There are two regions of integration, 
with characteristic scales of momenta p T ~ p M ~ m T and p T ~ m T , p^ ~ m M . In 
the case of the first region, which we can call the hard contribution, we can safely 
regard the muon mass and external momentum as small with respect to the integra- 
tion momenta, and expand the integrand in these small parameters. The resulting 
integral corresponds to a vacuum diagram shown in Fig. 3(a). In the second region 
we cannot neglect the external muon momentum, but now, since p M <C p T , we can 
expand the r loop propagators in Taylor series in p M , so that the integral factorizes 
into a product of two one-loop diagrams, shown in Fig. 3(b). This nice factorization 
of the relevant integration regions is only possible in regularization schemes, which 
do not introduce additional mass scales, such as dimensional regularization. 

After evaluation of these two sub-diagrams and renormalization of the electric charge 
of the muon, we find a finite result. Several of the terms we calculated agree with 
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Figure 3: Graphic representation of characteristic momentum scales in the two-loop 
diagram with a r loop. Solid and dashed lines denote, respectively, massive and massless 
propagators. 



formula (12) in [9]: 

. , I 2 / 4 ln/ 9 l4 131 l6 4/ 6 , , 

a 2 (l = m,./m T ) = 1 1 1 P H In / 

21 M/ T) 45 70 19600 99225 315 

~ (8n 3 + 28n 2 - 45)Z 2n+2 " n/ 2n+2 

- t 3 [(^ + 3)(2n + 3)(2n + 5)] 2+ n ^ (n + 3)(2n + 3)(2n + 5) ' 

The other two-loop diagram that we have to consider is the electron loop insertion 
in the photon propagator (Fig. 2). Here the situation is somewhat more involved, 
since there are now three integration regions. Introducing the obvious notation for 



~ m 



the integration momenta, p e ^, we have: (p e ~ p^ ~ m^), (p e ~ m e and p M 
and (p e ~ p^ ~ m e ). The first two regions correspond to known cases of the 
"large momentum expansion" [16]. They are depicted in Fig. 4(a,b). Figure 4(a) 
denotes a simple Taylor expansion of the electron propagators in m e , justified if both 
integration momenta are large. In Fig. 4(b) we have illustrated one of the two cases, 
where one of the electron lines cannot be expanded in m e , but where the integration 
factorizes and we have a product of one-loop diagrams. More exotic is the third 
case, Fig. 4(c), where the only scale of integration is m e , since all dependence on 
the muon mass and external momentum factorizes. Sub-diagrams of this type have 
been encountered in a different context in eikonal expansions [17, 18]. Integrals that 
arise in the present case are somewhat different and we give here some details of 
their evaluation. 

Since now the integration momentum p^ is much smaller than the external muon 
momentum p, we can expand muon propagators in p 2 . As a result we have to 
compute integrals of the form 

J( x J_ f d D Pfl d D p e 

{ai ' a2> W (p 2 )-(2p M .p)-(p 2 + m 2 )[(p M + p e ) 2 + m 2 ]- 



5 




(a) (b) (c) 

Figure 4: Integration regions contributing to the two-loop diagram with an electron 
loop, (a): Taylor expansion in m e . (b): region of soft electron loop momentum, (c): 
region of both momenta soft (double line denotes muon propagator expanded in the 
square of its virtual momentum). 



First, we use Feynman parameters to combine the last two propagators and integrate 
over p e . We get 

d D P , 



J(a 1 ,a 2 ) = -J^T(e)x e (l - x) 6 J 



(8) 



with m 2 = m\jx{\ — x). Next, using again Feynman parameters, we combine the 
first and the last term in the denominator. Finally, using 



I foo 

dX 



a«bP B(a,P)Jo [a + 6A] a+/3 ' 
and integrating over p M , we get 

r (f ) r (oi + f - 2 + 2e) r (2 - ai - f - e) r 2 (-1 + ai + f + e) 



(9) 



J(a 1 ,a 2 ) = 



2r(a 2 )r (2 - f - e) T(-2 + 2ai + a 2 + 2e) 



x- 



m 



4— 2ai— a2-4e 



m a 2 



7T 2 r ( ai + f - 2 + 2e) m^ai-aa-^ 



2 2a i+ 2a 2-3+2e r (a2±l) p ^ + + ^ p ^ _ |L _ e ) S i n ^ ( Ql + f + e ) m £ 2 ' 

We should mention that the integrand in eq. (7) could also contain products p ■ p e . 
It is possible to replace them by combinations of products of p^ ■ p e and p M • p using 
traceless products (see e.g. [17]). However, in the present case there are at most two 
powers of p e in the numerator and we can use the following simple formulas: 



ip-Pef 



P 2 PI 
D 



+ 



(P,-Pf 



'" n D(p e -p^f -p 2 e pl 



PP\ 
D 



{d-i){pD 



P-Pe 



pI 
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Adding the contributions of the three integration regions, we find, after renormal- 
ization, that the terms we obtained agree with formula (11) of [9]: 



ai(k = m e /m^) 



25 



71 



+ —k - - \nk + (3 + 4\nk)k 2 - -rr 2 k 3 
36 4 3 4 



7T 



+ 

+ E 



44 14 



T + f-f^k + 21n 2 k 



k 4 + —k 6 Ink 
15 



109 



n=2 



2(n + 3) 



n(2n+ l)(2n + 3) 



In k — 



225 

8n 3 + 44n 2 + 48n + 9 
n 2 (2n + l) 2 (2n + 3) 2 



k 2n+A 



3 Three-loop diagram with e and r loop inser- 
tions 



We now proceed to the actual focus of our work, the contribution of the three-loop 
diagram with e and r loop insertions, shown in Fig. 1. So far its contribution to a M 
has been evaluated only numerically [12]. 

Using the techniques described above, we can easily obtain an expansion of the 
Pauli part of this diagram with arbitrary accuracy. There are now three integration 
momenta and we have to consider five integration regions, combinations of the con- 
ditions described in the context of two-loop diagrams. Using the notation p T ,^,e for 
the integration momenta in the three loops, the regions we have to consider are: 

• Pr ~ ~ Pe ~ TTlr, 

• p T ~ m T ; p M ~ p e ~ m M , 

• p T ~ m T ; p M ~ m M ; p e ~ m e , 

• p T ~ m T ; p M ~ p e ~ m e , 



Calculations in each of these regions are analogous to the cases described in the 
previous section. For present purposes it is more than sufficient to retain the first 
three terms in the m 2 /m 2 expansion and two terms in m 2 /m 2 . After renormalization 
we find 



Cr pol -(e,r)^4f— ln^-— ) 
mi V135 m e 135/ 



mt f 229213 vr 2 37 , m T 1 , m T , m T m„ 3 , m„ 

+— £ + In — In — In — ^ + In — 

r< ^ 12348000 630 11025 m M 105 m M m 2 4900 m e 
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ml ( 1102961 4tt 2 398 m. 
mf ^~ 75014100 + 2835 ~ 297675 n m 
2 ml 4ir 2 ml 

j c c 

15 ml 45 m^m^ 
= 0.0005276(2), 



8 m r m T m u 524 m w \ 
m In — In — - 

945 m M m 2 e 297675 m e / 

(10) 



in agreement with the numerical evaluation. The error in the result is due to the r- 
lepton mass uncertainty. The leading-logarithmic term of this expansion corresponds 
to simply replacing a(q 2 = 0) by ot(mfy in the two-loop diagram with a r loop. We 
have included the last term, with odd powers of m e and m^, even though it is not 
relevant numerically. It illustrates typical contributions of the eikonal expansion, 
the only source of terms non-analytical in masses squared. 

We have checked eq. (10) by comparing it with an analytical integration, using the 
kernel function given in [13]. Terms quadratic in masses m TjjUie are in complete 
agreement. The odd powers of m e have not been included in [13]. 

With formula (10) the complete QED contribution to is now known analytically. 
In this particular case this result does not noticeably improve the accuracy of the 
QED prediction, since the error in (10) comes from the r lepton mass measurement. 
However, the technique presented here might facilitate other calculations. We have 
seen that a combination of large momentum, heavy mass, and eikonal expansions 
eliminates the need for numerical calculations and enables us to construct arbitrarily 
accurate expansions without knowledge of the exact result. 
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